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1.  Introduction 

The  main  objective  of  our  STTR  Phase  II  program  is  to  design  and  build  a  chip-scale  cost 
effective  reconfigurable  optical  device  that  relies  on  lasing  bistability  in  a  system  of  single  and 
coupled  microresonators  with  embedded  quantum  dots  (QDs). 

Relevance:  Military  decision  making  at  both  tactical  and  strategic  levels  crucially  depends  on 
capabilities  for  real  time  processing  of  large  amounts  of  information.  All-optical  platforms  are 
viewed  as  the  main  technology  capable  of  delivering  the  required  speed  for  information 
processing.  The  reconfigurable  photonic  device  being  developed  by  Hybrid  Photonics  will  play 
the  role  of  a  transistor  in  photonic  integrated  circuits  providing  all-optical  signal  processing  at 
speeds  exceeding  Tbits/second  and  low  power  consumption. 

Approach:  Two  different  systems  were  pursued  for  realizing  reconfigurable  optical  systems  that 
show  all-optical  bistability.  The  first  system  consists  of  coupled  microresonators.  Here 
microdisks,  microring,  vertical  cavity  structures  and  coupled  distributed  feedback  structures  were 
investigated  to  obtain  bistable  lasing  modes.  The  second  system  consisted  of  a  passive  microring 
resonator  coupled  to  two  bus  waveguides  which  have  an  integrated  gain  medium.  The  active 
material  in  all  of  these  systems  is  colloidal  quantum  dots.  While  the  present  demonstrations  were 
carried  out  using  CdSe/ZnS  dots  that  emit  in  the  red,  we  can  easily  extend  the  present  work  into 
the  infrared  using  PbSe  QDs. 

Summary:  The  biggest  challenge  in  the  phase  II  program  was  the  incorporation  of  colloidal  QDs 
into  host  materials  that  do  not  affect  their  optical  properties  while  being  able  to  be  patterned  into 
photonic  structures.  The  following  are  the  significant  outcomes  of  the  program: 

•  Successful  incorporation  of  CdSe  QDs  into  polymer  and  dielectric  host  and  realization  of 
devices  such  as  active  waveguides,  microdisk  and  microring  resonators,  surface  emitting 
microcavity  structures,  and  distributed  feedback  (DFB)  grating  structures 

•  Demonstration  of  gain  and  amplified  spontaneous  emission  (ASE)  from  a  vertical  cavity 
surface  emitting  structure. 

•  Demonstration  of  gain  and  lasing  like  emission  from  planar  DFB  grating  structures 
fabricated  using  softlithography. 

•  Development  of  fabrication  protocol  for  developing  surface  emitting  and  DFB  lasers 
using  colloidal  QDs  as  the  active  medium. 

•  Development  of  theoretical  models  for  coupled/single  resonator  bistable  lasing  structures. 
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2.  Coupled  Microresonator  System 

2.1  Theory  and  Software  Development 

Theory/software  component  under  this  approach  involved  the  development  of  simulation  tools 
for  modeling  the  lasing  characteristics  of  coupled  microdisks/ring  resonators.  We  also  performed 
finite  difference  time  domain  simulations  to  model  the  electromagnetic  field  distributions  in  the 
coupled  resonator  system. 

2.1.1  Modes  of  passive  disks 

In  the  course  of  the  project  we  developed  a  model  and  corresponding  software  for  computing 
modes  of  the  coupled  disks  structures  using  generalization  of  the  “constant-flux-mode”  approach 
[1,2],  which  was  applied  to  a  single  disk  problem  during  Phase  I  of  the  project.  We  consider  i  - 
th  disk  of  radius  Rt  and  refractive  index  nd  whose  center  is  positioned  at  r  =  r;  and  introduce 
internal  field  inside  the  i  -th  disk  as  a  linear  combination  of  Bessel  functions  of  the  first  kind, 


F(i)  =  V d(i)eim<pJ  (n,k  |r-r.|)  (1) 

int  /  j  m  m\d  ms  \  i  \ )  V  / 

m 

where  kms  is  the  eigenfrequency  of  the  coupled  system  (we  use  the  system  of  units,  in  which 

speed  of  light  in  vacuum  is  unity)  characterized  by  the  azimuthal  number  m  and  radial  number 
5  .  The  field  scattered  by  the  i  -th  disk  is  introduced  via  linear  combination  of  Hankel  functions 
of  the  first  kind,  Hm  (x) 


F„='Lby~''H.(r,ak\r-,\)  (2) 

m 

where  n0  is  the  refractive  index  of  the  external  medium,  and  A:  is  a  free  spectral  parameter 
representing  frequency  of  out  outgoing  stationary  lasing  output.  The  distinction  between  outside 
spectral  parameter  and  the  eigenvalues  of  the  system  (which  become  functions  of  £)  is  the 
characteristic  feature  of  the  constant-flux-mode  approach  [1,2].  Function  F  in  Eq.  (2)  represents 
normal  to  the  disk  component  of  the  magnetic  field  for  TE  polarization  and  that  of  electric  field 
for  TM  polarization.  The  field  incident  on  this  disk  consists  of  the  fields  scattered  by  all  other 
disks: 


j^i  m 

where  polar  angles  <pj  and  radial  coordinates  are  defined  in  local  coordinate  systems  associated 

with  each  j  -  th  disk.  In  order  to  be  able  to  use  boundary  conditions  at  the  rim  of  the  i  -  th  disk 

we  need  to  rewrite  the  scattered  field  in  the  coordinate  system  centered  at  the  i  -  th  disk.  This  is 
achieved  with  the  help  of  Graffs  formula  for  the  Hankel  function  [3] 
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r-r. 


n= ou 

)=Z". 


. |r  -  r  I). 


m(pi 


(4) 


where  Rjt  and  6fi  are  radial  and  polar  coordinates  of  the  /  -  th  disk  in  the  coordinate  system 

centered  at  j -th  disk.  Using  standard  Maxwell  boundary  conditions  at  the  rim  of  the  disk  we 
derive  for  the  expansion  coefficients  of  the  internal  field 


d(i)  =d(i)(k,k  )YY  KU)dU)H  ( nM)e 

m  m  V  ’  m,s  /  /  j  /  j  n  n  m-n  \  0  J1/ 


i(n-m)6ji 


(5) 


j^i  n 


where 


^m(kAm,s)  = 


2  in 


d,i 


™A  n0kmsHm  (kn0Rt )  J'm  (kmsndRi ) - kndJm  (kmsndRi ) H'm  {kn{)Ri ) 
2  i  1 


,  TE 


(6) 


”Ri  nAnsH,„  (MA  )  J'm  (KsndRi  )  -  knoJm  ( KSndRi  )  H'm  (MA  ) 


,TM 


and 


iA^kn"J'm  (ndkm,R,  )  -  ndK,Jl,  (ndkm,R,  )Jm  (hlA  )] 

i!I!A[knaJ’m  (kn„R, )  J,  (nXA )  -  n,kmJ',  ("A,*,  K  {k"A, )] 

a 


id) 


In  the  case  of  only  two  disks  we  can  choose  the  axis  of  the  coordinate  system  along  their  centers 
so  that  6n  =O;02l-7r  and  Eq.  (5)  takes  the  following  form 


<C  =  oi"  A  k,  )£  *-f  (n,kRn ) 

n 

421  =  «i21  (A K  >E  (n„kRn ) 


For  high-Q  resonances,  the  largest  contribution  to  the  sum  over  modes  in  Eq.  (8)  comes  from  the 
term  with  n  =  -m  .  In  this  resonant  approximation  we  have 


dm  =  d(x\k,k  )/c(2)L  d(2); d(2)  =d(2\k,k  )jc(1)L  d(l) 

m  m  V  ’  s  ’  m  2m  -m  ?  -m  m  V  ’  s'  m  2m  m 

d(l)  =a(l\k,k)K(1\  d(2)-d{2)=a(2\k,k)K{%  dil) 


(9) 


where  we  introduced  a  coupling  constant 

^2 m  =  H 2m  (n0kR\2  ) 


(10) 
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The  system  of  equations  (9)  describes  two  pairs  of  degenerate  modes,  whose  eigenfrequencies 
are  given  by  equation: 


1  =  a £1}  (k,  ks  (k,  ks )/ci1)42)  (t2m  )2 


(11) 


This  equation  has  two  solutions  in  the  vicinity  of  m-th  mode,  which  we  will  designate  as 
kmA  and  km  l .  In  order  to  look  for  these  solutions  numerically  we  exploit  the  same  algorithm  that 

we  used  for  a  single  disk  problem  in  the  Phase  I.  In  order  to  find  the  starting  points  we  can  try  to 
use  the  following  approach.  If  the  disks  are  identical  (or  almost  identical),  we  present  the 
scattering  coefficients  in  the  vicinity  of  the  single  disk  solutions  as 


«i"(/v,/yK 


.(1)  _ 


/  m,s 


k~ks  +  iy, 


(1) 


g(M>:21 = .  ,2, 

k-k+  iyl  , 

s  /  m,s 


Eq.  (11)  then  takes  the  form  of 


/  m,s 


-ir(2) 

/  m,s 


k-K+ifml  k-K+iy, 


(t,J= 1 


[*-*,+  irSs  ][k-k,+  irZ  ]  =  -rX  (t,m  f 


(12) 


(13) 


which  can  be  solved  to  give 


-iym  -  k{1)  +  iy(2)  )2  -  4y(1)  y(2)  ( t2  )2 

/  m,s  m,s  /  m,s  J  /  m,s/  m,s  \  2m  ) 

-iym  -k(2)  +iy(2)  )2  -4r(1)  r(2)  (t7  )2 

/  m,s  m,s  /  m,s  J  /  m,sf  m,s  \  2m  / 


(14) 


These  values  are  used  as  initial  approximations  for  each  of  the  two  possible  solutions.  Results  of 
numerical  solution  of  Eq.  (11)  are  shown  in  Fig.  1 
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Fig.  1.  Dependence  of  the  real  (left)  and  imaginary  (right)  parts  of  two  modes  of  the  two-coupled  disks  system  on  the 
distance  between  the  disks.  Radii,  R ,  of  the  both  disks  are  the  same. 

Taking  into  account  the  double-degeneracy  of  the  modes  of  the  two-disk  system,  one  can  define 
four  linearly  independent  modes  as  follows: 


(1)  :  d(l\  =  N  ,d(l)y2),L  ;  d(2) .  =  N  . ;  d(l) .  =  0;  d(2),  =  0; 

V  /  m,  1  m,l  m,l  m,\  2m  ’  -m,  1  m,  1’  -m,  1  ’  m,  1  ’ 

(2):d(l\=N  id(1\/c(2)0L  ;d(2\=N  d(l) .  =  0;  d(2\  =  0; 

\  f  m, 2  m,2  m,2  m,2  2m  ’  -m,2  m,2  ’  -m,2  ’  m,2  ’ 

(3)  :  =  Nmla^>t2.;  <>  =  =  0;  =  0; 

(4)  :  <7™.  4  =  NmAa'"2K'»2t2m,d™  =  JV„  4;<>  =  0;  </«>„  =  0; 

where  we  introduced  an  abbreviated  notation 


«(U)=«(1’2) 
m,i  m 


(k’km,i)’ 


*-0.2)  =  *-0,2) 
m,i  m 


(15) 

(16) 

(17) 

(18) 

(19) 


The  fields  in  the  respective  modes  are 
Mode  1 


Mode  2 


F(l)=R~lN  a(l)K(2)H  (nkR  )e‘mtpJ  ink  Ir-rl) 

1  1,1  Jv  ^  m^m^mF12m\n0K1\2  )K  |*  *l|/ 

*u  . (»A,  |r-r,|) 

*u  («A,.  -r2|) 


(20) 


(21) 


Mode  3 
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N’  =  iX.  {ndkm2 

F$=R-'K„2e‘-tJ_m(ndK3\r-r^ 


r-r 


i) 


(22) 


Mode  4 


FI  =  R-'X,jS 1VS^2.  («.**,,  K""V  .  |r  -  r,  |) 

^,=«-,^^.(»A.2|r-r2|) 


(23) 


where  i?is  characteristic  size  of  the  resonators  used  in  the  definition  of  the  normalization 
constants  Nm  i  (see  Eq.  (28)  below).  If  disks  are  identical,  it  is  their  radius.  Now,  we  need  to 

define  the  inner  product  of  the  modes.  To  this  end  we  first  introduce  two-dimensional  vectors 


and  their  adjoints 


f(2) 

\ri  J 


(24) 


(25) 


where  : 


'F?\r,(p)  =  F?\r,-cp) 


(26) 


Now  the  inner  product  is  defined  as 

2k  ^2  2.K 

<N|F  )  =  {  |  fF7(1) (r, (p)F^ (r, (p)rdrd(p  +  J  J  fFy(2) (r, <p)F^ ( r , (p)rdrd(p  (27) 

0  0  0  0 


One  can  prove  that  the  modes  introduced  in  Eq.  (20)  through  Eq.  (23)  are  orthogonal  with 
respect  to  the  inner  product  defined  by  Eq.  (27).  Normalization  of  the  modes  is  defined  by  the 
following  normalization  condition 


J  K“i2 )]’ 

1  2  1  2  1 

X»i 2\\/m(ndKAp)\  PdP  +  Vl\\_JAndkm,R2p)\  Pdp\  =  1, 
0  0  J 


(28) 


where  ul-Rl/R;u2  =  R2/R 
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2.1.2  Lasing  equations  and  analysis  of  single  mode  operation 

The  relatively  low  Q-factor  of  the  coupled  disk  system  and  the  presence  of  degenerate  modes 
make  standard  lasing  rate  equations  unsuitable  for  modeling  the  structure  under  consideration. 
Therefore,  parts  of  our  efforts  were 
devoted  to  derivation  of  the  system  of 
equations  for  amplitudes  of  lasing 
modes.  For  details  on  the  derivation  of 
the  lasing  equations,  please  see  the 
appendix.  These  equations  were  used  to 
analyze  the  single  mode  operation  of  the 
double  disk  laser  as  well  as  compute 
coefficients  of  nonlinear  inter-mode 
coupling,  which  are  essential  for 
determination  of  the  system  parameters 
suitable  for  the  bistable  behavior.  The 
results  of  our  modeling  are  shown  in  Fig.  2.  The  lasing  frequencies  of  the  two  non-degenerate 
modes  as  functions  of  the  frequency  of  quantum  dot  (QD)  transition  in  the  regime  when  only  one 
of  the  modes  is  assumed  to  attain  lasing  is  shown  in  Fig.  2(a).  Figure  2(b)  presents  dependence 
of  the  single  mode  lasing  threshold  for  the  same  two  modes  also  as  functions  of  QD  frequency. 

2.1.3  Finite  Difference  Time  Domain  Simulation  of  Coupled  Microring  Resonator 

We  performed  finite  difference  time  domain  simulations  of  the  coupled  microrings  and 
established  the  mode  splitting  that  happens  in  these  resonators.  Transmission  spectra  at  the  two 
ports  are  shown  in  Fig.  3  clearly  showing  the  mode  splitting.  Also  shown  in  Fig.  3  is  the  electric 
field  distribution  in  the  coupled  ring  resonator  system  at  the  resonances. 


Atomic  frequency 


Fig.  2  Dependence  of  the  single  mode  lasing  frequency  and  threshold  for 
two  modes,  as  functions  of  the  QD  frequency  for  two  inter-disk  distances, 
d ,  measured  in  disk’s  diameter. 


Fig.  3.  Transmission  characteristics  of  the  coupled  microring  resonator  showing  the  splitting  of  modes.  Also  showin  on 
the  right  are  the  electric  field  profiles  in  the  coupled  resonators  at  the  three  spectral  positions  (a),  (b)  and  (c). 
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2.2  Fabrication  of  Active  Micro-resonators 


SiNx<*«l 


Si02 


Si 

Fig.  4.  Schematic  drawing  of  the  active 
microdisks  realizing  using  colloidal 
quantum  dots  embedded  in  a  SiNx  host. 


2.2.1  Dielectric  Microdisks  with  Quantum  Dots 

On  the  experimental  front,  we  fabricated  active  and 
passive  microdisks.  Although  we  had  originally  proposed 
to  use  SU8  as  the  host  medium  to  embed  the  colloidal 
quantum  dots  to  realize  the  active  disks,  the  QC  team 
recently  introduced  SiNx  as  a  host  matrix  with  a 
significant  advantage  of  having  very  small  detrimental 
effect  on  the  luminescence  properties  of  the  embedded 
dots.  Schematic  drawing  of  the  active  microdisk  is  shown 
in  Fig.  4.  The  SiNx  disk  is  fabricated  on  a  Si02  buffer 
layer  which  acts  as  the  low  index  bottom  cladding  layer  as 
well.  The  SiNx  disks  vary  in  diameter  (5  pm  to  20  pm)  and  are  0.5  pm  thick.  The  entire  structure 
is  realized  on  Si  substrate  and  grown  using  a  combination  of  plasma  enhanced  chemical  vapor 
deposition  and  spin  coating. 

Following  the  deposition  of  the  different  layers  in  the 
following  order:  Si02,  SiNx,  QDs,  SiNx,  the  microdisks  are 
patterned  using  electron  beam  lithography  using  PMMA  as  the 
electron  beam  resist.  The  patterned  PMMA  is  then  used  as  the 
etch  mask  to  etch  the  microdisks  onto  the  SiNx  layer  using 
reactive  ion  etching.  Scanning  electron  microscope  images  of 
coupled  passive  disks  and  single  active  disks  are  show  in  Fig.  5. 

CdSe  colloidal  QD  with  emission  maximum  at  600  nm 
wavelength  was  used  as  the  active  component  in  these  disks. 

These  QDs  were  deposited  using  spin  coating.  It  was  observed 
that  the  sidewalls  of  the  microdisks  fabricated  using  the 
above  technique  suffered  from  significant  sidewall 
roughness  which  in  turn  affected  the  luminescence 
properties. 

Focused  ion  beam  etching  was  explored  as  an 
alternate  technique  to  realize  smooth  sidewalls  for  the 
microdisks.  This  was  carried  out  at  Brookhaven  National 
Labs  -  Center  for  Functional  Nanomaterials  (Fig.  6). 


Fig.  5.  SEM  image  of  active  microdisk 
embedded  with  red  emitting  CdSe 
colloidal  quantum  dots  (left)  and  coupled 
microdisks  with  500  nm  spacing  (right). 


Fig.  6.  SEM  image  of  active  microdisks 
patterned  using  focused  ion  beam  etching.  As 
can  been  seen,  the  sidewalls  are  far  soother 
than  disks  patterned  using  RIE, 


2.2.2  Polymer  microdisk/rings 

In  addition  to  the  all¬ 
dielectric  microdisks,  we 
have  also  investigated 
single  and  coupled 
microdisk/ring  structures 
embedded  with  colloidal 
CdSe  quantum  dots  and 
Rhodamine  6G  organic 
dye.  The  choice  of  this 
dye  was  motivated  by 
the  high  photostability 


Fig.  7.  (a)  SEM  image  of  polymer  microdisks  embedded  with  Rhodamine  6G.  (b)  Sidewall 
images  of  a  coupled  disk  system  and  (c)  Photoluminescence  spectra  collected  from  the 
microdisk  resonator. 
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and  high  quantum  yield  (0.95).  This  approach  provides  an  alternate  approach  for  realizing  low- 
cost  microlasers.  The  high  quantum  yield  of  the  dye  as  compared  to  the  colloidal  QDs  was  the 
primary  motivation  for  pursuing  this  approach.  We  expect  the  dye  doped  material  to  exhibit 
lasing  at  lower  thresholds  than  the  quantum  dot  system.  The  dye  was  dissolved  along  with  SU8 
(Microchem  Inc.)  and  thinned  using  a  thinner  to  obtain  one  micron  thick  films.  The  dye  doped 
films  were  later  patterned  using  either  EBL  or  photolithography.  Scanning  electron  microscope 
images  of  the  dye  doped  disks  fabricated  in  SU8  is  shown  in  Fig.  7  along  with  the  emission 
obtained.  Once  again  we  see  the  surface  roughness  of  the  microdisks  to  be  a  hindrance  in 
observing  the  WGMs.  It  should  however  be  noted  that  in  both  the  dye  doped  samples  as  well  as 
in  the  QD  embedded  disk  samples,  light  was  collected  in  the  far-field  and  this  often  deteriorates 
the  signal  integrity.  A  side  out-coupling  geometry  is  currently  being  designed  to  enhance  the 
coupling. 

We  fabricated  Rhodamine  6G  infiltrated  polymer 
microrings  using  the  same  fabrication  protocol  established 
for  the  microdisk  resonators.  Scanning  electron  microscope 
images  of  single  and  coupled  microring  resonators 
fabricated  using  EBL  are  shown  in  Fig.  8.  Far-field 
emission  measurements  showed  very  weak  PL  signal  since 
the  amount  of  dye  in  the  rings  is  far  smaller  than  that  in 
the  disks.  As  a  natural  direction  in  this  project,  we  are  now 
developing  structures  with  waveguides  coupled  to  the 
rings. 


. . -  IM  ■ 

Fig.  8.  SEM  images  of  single  and  coupled 
microring  resonators  infiltrated  with 
Rhodamine  6G. 


A  vertical  coupling  geometry  was  used  to  fabricate  the  passive 
waveguides  coupled  to  the  active  microrings.  Schematic  of  the  flow 
diagram  for  fabricating  the  structure  is  shown  in  Fig.  9.  The  passive 
waveguides  were  fabricated  at  QC-CUNY  using  electron  beam 
lithography.  This  was  followed  by  the  deposition  of  a  top  cladding  layer 
of  silicon  dioxide,  the  thickness  of  which  determines  the  coupling 
between  the  passive  waveguides  and  the  ring.  Finally  active  microrings 
were  defined  on  our  standard  silicon  nitride  structure  embedded  with 
QDs  or  dye  using  FIB  technique  at  Brookhaven.  The  major  fabrication 
issue  we  had  in  this  context  was  the  precise  alignment  of  the  passive 
and  active  waveguides. 
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Fig.  9.  Schematic  flow  diagram  for 
the  fabrication  of  the  vertically 
coupled  microring  system. 


2.2.3  Vertical  coupled  microcavities 

Recently  Menon’s  group  at  QC-CUNY  has  developed  a  fabrication  protocol 
that  uses  plasma  enhanced  chemical  vapor  deposition  (PECVD)  compatible 
with  colloidal  quantum  dots  without  deteriorating  the  optical  properties  of 
the  dots.  In  the  present  program  that  approach  has  been  utilized  to  develop 
coupled  vertical  microcavities  to  realize  bistable  lasing.  The  entire  cavity 
structure  is  fabricated  using  PECVD  and  the  quantum  dot  layers  are 
deposited  via  spin  coating  or  dip-coating.  Shown  in  Fig.  10  is  the  schematic 
of  the  coupled  vertical  microcavity  structure.  The  number  of  periods  of  the 
middle  DBR  was  controlled  to  manipulate  the  coupling  between  the 
microcavities. 


Fig.  10.  Schematic  of  a 
coupled  microcavity. 
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Fabrication 

Using  a  combination  of  nitrous  oxide,  silane,  and  ammonia,  we  deposit  layers  of  silicon  oxide 
(Si02)  and  silicon  nitride  (S^N^.  Plasma  enhanced  chemical  vapor  deposition  is  used  for  the 
growth  of  these  dielectric  layers.  The  SiC>2  and  Si3N4  layers  are  estimated  to  have  refractive 
indices  of  1.45  and  1.77,  respectively.  A  low  temperature  growth  technique  was  developed  to 
deposit  the  dielectric  layers  on  top  of  the  active  QDs  to  avoid  the  QDs  from  getting  affected  due 
to  the  growth  process.  This  was  an  important  process  development  in  the  program  since  it 
established  the  protocol  to  introduce  active  material  between  dielectric  layers  without  having 
detrimental  effect  on  the  optical  properties  of  the  QDs. 


Characterization  and  Modeling 

By  controlling  the  number  of  periods  of  the  middle 
the  coupled  microcavities.  Additionally,  the 
reflectivity  and  transmission  spectrum  shows  the 
split  modes  arising  due  to  the  interaction 
between  the  microcavities  as  seen  in  Fig.  11. 
The  top  panel  of  Fig.  11  shows  the  split  modes 
arising  when  the  number  of  intermediate  DBR 
periods  is  small  (3.5)  and  the  bottom  panel 
shows  the  reflectivity  spectrum  showing  lesser 
splitting  of  the  modes  since  the  number  of 
intermediate  DBR  periods  is  larger  (10.5).  The 
former  case  corresponds  to  strong  interaction 
between  the  two  cavities.  The  experimentally 
obtained  reflectivity  data  is  compared  to 
simulations  based  n  transfer  matrix  method.  The 
corresponding  electromagnetic  field  profiles  are 
shown  in  the  right  panel  of  Fig.  11.  As  can  be 
seen,  it  is  possible  to  precisely  control  the  field 
distribution  in  these  coupled  resonator  systems 
by  engineering  their  coupling  strengths.  The 
electromagnetic  field  simulations  were  also 
carried  out  using  transfer  matrix  method. 
Following  this  we  carried  out  time  resolved  and 
steady  state  photoluminescence  measurements 
on  the  emission  characteristics  of  the  QDs 
embedded  in  the  coupled  microcavity  structures 
(Fig.  12).  A  Purcell  enhancement  of  2.5  was 
observed  for  the  system  of  coupled  cavities 
embedded  with  CdSe  QDs  in  one  of  the  cavities. 

3.  Single  Resonator  System 


DBR,  one  can  control  the  field  distribution  in 


Fig.  11.  Coupled  microcavities  reflectivity  with  fitted  spectra, 
side-by  with  electric  field  intensity  profile.  Separations  between 
the  two  cavities  are  3.5,  and  10.5  periods  for  top,  and  bottom 
graphs  respectively.  The  field  intensity  is  represented  by  the 
color,  vertical  and  horizontal  axis,  are  position  within  the 
structure  and  wavelength,  respectively. 


Fig.  12.  Photoluminescence  (PL)  from  QDs  in  coupled 
microcavity.  (Left)  Steady  state  and  (Right)  Time  resolved  PL. 


3.1  Theory  of  linear  bistability  mechanism 

Significant  efforts  have  been  devoted  this  year  to  exploring  a  new  type  of  optic  bistability  in 
lasers,  which  does  not  rely  on  mode  competition  or  saturable  absorbers,  and  is  expected, 
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therefore,  to  show  bistable  behavior  at  lower  power 
consumptions  and  with  faster  switching  times.  This 
mechanism  of  bistability  is  realized  in  structures  shown  in 
Fig.  13,  where  two  active  waveguides  (length  of  active  portion/,) 
coupled  to  each  other  via  a  disk  or  ring  resonator  (circumference  Lr , 

index  nr ).  The  feedback  necessary  for  lasing  is  provided  in  this 

by  combination  of  two  factors:  mirrors  at  one  end  of  each 
waveguide  and  coupling  of  the  waveguides  via  the 
resonator,  which  redirect  optical  signal  reflected  from  one 
waveguide  mirrors  to  another  waveguide  closing  the 
loop.  Normal  modes  of  this  structure  are  defined  by  the 
following  equations  for  the  real  parts  of  the  propagating  constants  in  the  waveguides,  q 


are 

refractive 

structure 


Fig.  13.  A  schematic  of  the  structure  suggested 
for  realization  of  the  linear  bistability 
mechanism 
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and  their  imaginary  parts  Sj'/2’ 


Eq.  (29)  indicates  that  the  disk  does  not  affect  the  modal  structure  of  the  system,  which  is  similar 
to  that  of  an  effective  Fabry-Perot  resonator  of  length  2 L.  Its  role  is  to  determine  the  modal 
losses,  characterized  by  S„’2) ,  via  resonator  form-factor  S(oj)  defined  as 


S  =  - 
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a  r  \k\ 


1  +  a4 1/?|4  -  2a2 1/?|2  cos 


o  0J 

2n — 


co. 


r  J 


where  cor  =  {2k c)  /  {nrLr )  ( c  is  speed  of  light  in  vacuum),  a 2  describes  round-trip  losses  in  the 

resonator,  p  and  k  are  coupling  parameters  obeying  |/j>|2  +  |/c|2  =  1 .  The  resonator  also  introduces 
a  phase  shift  between  the  fields  in  the  waveguides  defined  by 


V 


0/2- tan  1 


|jo|2  a  sin® 
l-|jo|2  a  cos® 


where  O  =  n  coL  /  c . 

r  r 

We  confirmed  these  conclusions  about  the  modes  of  our 
structure  by  direct  numerical  simulation  of  the  waveguide-resonator 
system  using  COMSOL  simulation  platform.  The  results  of  the 
simulations  showing  the  structure  of  the  modes  is  shown  in  Fig.  14. 

Using  found  modes  we  developed  a  theory  of  a  single-mode  lasing  in 


#  % 
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Fig.  14  Field  distribution  in  a 
coupled  waveguide-resonator 
structure  in  one  of  its  modes. 
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this  structure.  According  to  the  theory,  the  stationary  lasing  frequency,  Q  and  respective 
intensity  I  are  defined  by  the  system  of  equations 


(31) 


where  v  is  the  frequency  of  the  “cold”  cavity  mode,  a>0  and  y±  are  the  frequency  and  the 
linewidth  of  the  lasing  atomic  transition,  g  (Q)  is  the  gain  Lorentzian,  P  is  pumping  rate,  and 
Fl2(I,Q.)  are  integrals  of  the  modal  functions  U  (z) 


(32) 


taking  into  account  gain  saturation  through  saturation  parameter  R  =  /|t/(z)|2g(Q)(Ais  mode 


normalization  parameter).  We  established  that  under  certain  conditions  controllable  by  the  design 
of  the  structure,  Eqs.  (3.1.2)  have  multiple 


power  Px  <  P2  <  P3  Below  the  lowest  one  no  lasing  £ 
takes  place,  for  PX<P  <P2 one  has  a  standard  2 
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single-mode  laser.  However,  when  P2<P<P3 the  j 


non-lasing  solution  becomes  one  of  two  stable  "  o.o  L  - '  -  — . . . , 
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solutions  so  that  depending  on  initial  conditions  the  P  pumpjng 

system  can  either  laze  or  not.  This  regime  is  ideal  Fig.  15.  (a)  Dependence  of  stationary  lasing  frequency  and 

for  developing  an  optical  logic  element,  in  which  intensity  on  pump  power  and  (b)  Temporal  evolution  of 
.  .  lasing  frequency  and  intensity 

lasing  regime  would  correspond  to  1  and  non¬ 
lazing  to  zero.  The  switching  between  the  regimes  can  be  realized  using  injection  locking  of  the 
initial  frequency.  Further  increase  of  the  pumping  results  in  the  zero  intensity  solution  losing  its 
stability.  However,  in  its  place  another  lasing  solution  with  a  frequency  and  intensity  different 
from  the  first  one  becomes  stable,  so  that  now  one  has  two  stable  lasing  solutions  characterized 
by  different  frequencies  and  intensities.  The  bistability  discussed  here  is  not  characterized  by  a 
typical  hysteresis-like  dependence  of  intensity  versus  pumping,  which  is  typical  for  other  types 
of  bistability. 
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Fig.  16.  Top  panel  shows  SEM  images  of  active-passive  integrated  polymer 
waveguides.  The  active  waveguides  consist  of  CdSe  quantum  dots  dispersed 
in  SU8.  Bottom  panel  shows  CCD  images  of  optical  pumping  experiment 
clearly  demonstrating  waveguiding  and  minimal  scattering  from  due  to 
clustering  of  Quantum  dots. 


3.2  Fabrication  of  active/passive  integrated  waveguides 

We  fabricated  colloidal  QD 
embedded  waveguides  to  study  the 
transmission  losses.  The 
waveguides  were  fabricated  using 
CdSe  colloidal  quantum  dots 
dispersed  in  SU8.  The  patterning 
was  done  using  electron  beam 
lithography.  Figure  15  shows  SEM 
images  of  one  such  waveguide 
embedded  with  CdSe/ZnS 
core/shell  quantum  dots  emitting  at 
600  nm.  The  L-shaped  quantum 
dot  embedded  waveguide  was 
optically  excited  using  an  Ar-ion 
laser  (k=  488  nm)  from  one  of  the 
ends  and  the  emission  from  the 

quantum  dots  was  found  to  be  clearly  waveguiding.  Based  on  the  CCD  image  (Fig.  15  bottom 
panel)  it  can  be  concluded  that  the  scattering  induced  loss  is  very  minimal  as  we  do  not  see  much 
scattered  light  from  the  side  walls.  Most  of  the  light  emission  is  observed  from  the  ends 
indicating  good  waveguiding.  Furthermore,  the  waveguides  did  not  show  much  losses  at  the 
bends.  This  is  a  very  promising  result  towards  realizing  active  photonic  structures  using  colloidal 
quantum  dots  which  can  be  integrated  with  passive  structures. 

4.  Enhanced  Gain  from  QDs  in  Microcavity 

Using  a  combination  of  PECVD  and  dip  coating,  we  fabricated  one-dimensional  microcavity 
structures  as  shown  schematically  in  Fig.  17. 

The  dielectric  layers  of  the  MC  structure  were  fabricated  using  plasma  enhanced  chemical  vapor 
deposition  (PECVD).  The  MC  consists  of  two  distributed  Bragg  reflectors  (DBRs)  separated  by 
a  thin  layer  of  CdSe/ZnS  (core/shell)  QDs.  The  DBRs  consisted  of  alternating  layers  SiC>2  and 
ShN/t.  The  bottom  DBR  consisted  of  10  pairs  of 
ShNzt/SiCh  deposited  on  a  glass  substrate.  The  cavity 
layer  consisting  of  a  QD  film  was  then  deposited  by 
dip  coating  in  a  solution  of  QDs  in  toluene  of  10 
mg/ml  concentration,  and  pulled  at  a  rate  of  1 
mm/min.  The  resulting  film  was  estimated  to  be 
~200  nm  thick.  The  top  DBR  consisting  of  9  pairs 
was  then  deposited  on  top  of  the  QD  layer 
(schematic  shown  in  the  inset  of  Fig.  1).  The 
reflectivity  spectrum  of  the  resulting  MC  at  normal 
incidence  was  measured  using  a  tungsten-halogen 
light  source  coupled  to  a  CCD  based  spectrometer 
via  an  optical  fiber.  The  reflectivity  spectrum 
consists  of  a  stopband  with  a  transmission  window 
centered  at  630nm  and  a  full-width-half-maximum 
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Fig.  17  (Inset)  Schematic  of  the  MC  structure.  The 
normal  incidence  reflectivity  of  the  MC  is  plotted 
together  with  the  PL  spectrum  of  a  solution  of  QDs  in 
toluene. 
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(FWHM)  of  ~10  nm,  which  spectrally  overlaps  with  the  photoluminescence  (PL)  spectrum  of  the 
QDs,  as  shown  in  Fig.  17.  A  control  sample  consisting  of  a  QD  layer  film  between  two  SiC>2 
layers  was  used  for  comparison. 


Steady-state  PL  was  used 
optically  pumped  with  a 
continuous-wave  (CW) 
argon-ion  laser,  and  the 
spectra  were  collected 
using  a  fiber-coupled 
CCD  based  spectrometer. 
The  FWHM  of  the  PL 
narrows  to  10  nm 
compared  to  the  control 
sample  line  width  of  32 
nm,  as  shown  in  Fig.  18a. 

The  PL  peak  wavelength 
of  the  MC  is  slightly  red- 
shifted  compared  to  the 
control  sample  due  to  the 
cavity  resonance  being  at 
a  higher  wavelength  than 
then  bare  QDs  (see  Fig. 
17).  We  observe  the  PL 
is  cavity  linewidth 
limited  (10  nm)  and  peak 
is  spectrally  located  at 
the  cavity  resonance 
(-625  nm,  see  Fig.  17) 
indicating  the  QD 
emission  is  coupled  to 
the  cavity  mode.  To 
demonstrate 

enhancement  of  PL, 
optical  power  dependent 
measurements  were 

performed.  The  PL 
intensities  for  the  MC  is 
higher  than  the  control 
for  all  powers  measured, 
however,  no  threshold 
behavior  is  observed  (see 
Fig.  18b). 

Furthermore,  time- 

resolve  PL  was  carried 
out  to  determine  the 


to  determine  the  effects  of  the  cavity  on  the  QDs.  The  QDs  were 


Fig.  18.  (a)  Normalized  PL  for  the  MC  and  control  sample.  The  PL  of  the  MC  sample  is 
narrower  (10  nm  FWHM)  than  the  control  (32  nm  FWHM).  (b)  Power  dependent  PL  peak 
intensities  show  the  MC  structure  has  greater  intensity  than  the  control  sample,  (c)  Time- 
resolved  PL  measurements  are  plotted  for  the  control  (red)  and  MC  (black)  sample. 


Fig.  19..  (a)  PL  spectra  for  pump  powers  varying  from  0.1  mW  to  15.5  mW.  The 
biexciton  peak  begins  to  appear  with  2.0  mW  pump  power,  (b)  The  biexciton  peak 
power  is  plotted  against  the  time- averaged  input  power  to  show  the  threshold  behavior. 
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radiative  lifetime  of  the  QDs  confined  in  the  MC.  To  this  end,  a  time-correlated  single  photon 
counting  technique  using  a  pulsed  467  nm,  <  200  ps  pulsewidth  laser  was  used.  The  lifetime  of 
the  QDs  in  the  MC  was  found  to  be  8.9  ±0.3  ns,  while  the  control  sample  was  measured  to  be 
1 1.8  ±  0.6  ns.  The  MC  structure  yielded  a  Purcell  factor  ~1.3.  The  time-resolve  PL  is  shown  in 
Fig.  18c. 

Power  dependent  PL 
measurements  were 
repeated  at  higher 
intensities  using  the 
second  harmonic  (532 
nm)  of  a  Nd:YAG  laser. 

Amplified  spontaneous 
emission  (ASE)  of  the 
biexcitons  begin  to 
appear  at  a  power  of  2.0 
mW,  as  shown  in  Fig. 

19(a).  The  PL  intensity 
of  the  biexcitons  exhibit 
a  threshold  behavior 
which  is  shown  in  Fig. 

19(b).  The  angular 
dispersion  allows  for 
tuning  the  cavity  mode 
to  be  in  resonance  with 
the  ASE  peak  by 
changing  the  collection 
angle.  Furthermore, 

angle  dependent  PL  show  the  effects  of  angular  dispersion  of  the  MC  mode.  The  biexciton  PL  is 
in  resonance  with  the  cavity  at  a  5°  angle  to  the  normal.  The  enhancement  of  the  gain  of  the 
biexcitons  is  clearly  seen  from  the  increase  PL  intensity  at  resonance  (Fig.  20(b)).  This  effect 
yield  highly  direction  emission  when  compared  to  the  Lambertian  pattern  as  shown  in  20(a). 

This  work  demonstrates  the  enhancement  of  gain  in  an  all-dielectric  MC.  Time-correlated  single¬ 
photon  counting  technique  was  used  to  measure  increase  of  33%  in  the  radiative  rate  of  the  QDs 
inside  the  MC.  Spectroscopic  measurements  show  the  enhanced  PL  power  of  the  QDs  in  the  MC 
compared  to  the  control  sample.  Pumping  at  higher  powers  enhance  the  gain  of  the  biexcitons. 
The  enhancement  of  light-matter  interaction  is  such  system  can  be  further  developed  into  and 
have  implications  in  low  cost  optically  active  devices. 


Fig.  20.  (a)  PL  emission  pattern  of  biexciton  ASE  in  MC.  The  emission  is  highly 
directional  when  compared  to  the  Lambertian  emission  pattern,  (b)  The  gain  of  the 
biexcitons  is  enhanced  when  in  resonance  with  the  microcavity. 


5.  Gain  measurements  using  variable  stripe  length  measurements 

To  better  understand  the  gain  characteristics  of  the  CdSe  QDs  used  in  the  device  demonstrations 
discussed  above  we  carried  out  gain  measurements  using  the  variable  stripe  length  technique. 
Here,  the  NdYAG  laser  was  focused  down  to  the  sample  using  a  cylindrical  lens  and  the  stripe 
length  was  determined  using  a  beam  blocker  with  a  micrometer  attached  to  it  to  vary  the  stripe 
length.  The  QDs  were  deposited  using  dip  coasting  technique  similar  to  the  MC  structures 
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Fig.  21  Results  of  variable  stripe  length  measurements  carried  out  on  the  CdSe  ODs. 


Wavelength  (nm) 

Fig.  22  Steady  state  PL  measurements  showing  significant 
spectral  narrowing  inside  the  microcavity. 


Wavelenght  (nm) 


Fig.  23  Power  dependence  of  PL  emission.  The  inset 
shows  the  directionality  in  emission. 


discussed  above.  The  gain  spectrum  was  measured  for  different  stripe  lengths  and  the  result 
obtained  is  shown  in  Fig.  21.  The  maximum  gain  observed  was  ~  1.8  x  10~3  pm1. 

Similar  devices  were  also  fabricated  using  PbSe  QDs  emitting  in  the  1550  nm  wavelength  range. 
The  MC  was  once  again  fabricated  via  PECVD  and  the  QDs  were  deposited  via  spin  coating. 
Results  of  the  near  IR  emitting  microcavity  device  is  shown  in  Fig.  22.  It  is  clearly  seen  that  the 
emission  spectrum  narrows  considerably  when  the  QDs  are  embedded  in  the  MC.  Although  gain 
was  not  observed  in  these  QDs,  we  do  see  clear  evidence  of  enhancement  of  spontaneous 
emission  due  to  the  Purcell  effect  (Fig.  22).  The  power  dependent  spectrum  is  shown  in  Fig.  23. 
The  inset  of  Fig.  23  shows  the  directionality  in  emission  observed  from  the  PbSe  QDs  clearly 
indicating  coupling  of  the  emission  into  the  cavity  mode. 
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6.  Bistability  in  nonlinear  grating 


We  have  looked  into  the  possibility  of 
observing  bistability  using  a  grating 
structure  for  feedback.  A  schematic  of  the 
lasing  structure  based  on  a  distributed 
feedback  geometry  is  shown  in  Fig.  22. 

Calculations  are  carried  out  for  small-signal 
amplification.  In  the  small-signal  regime, 
the  optical  power  is  too  small  to  affect  the 
gain  value.  Amplifier  gain  G  -  Iout  jlm  is 

calculated  for  a  A/4  phase  shifted  DFB.  A 
schematic  of  the  such  an  experiment  is 
depicted  in  figure  1 . 

The  parameters  used  were  kL  =  0,  aintL  -  0 . 


Fig.  22.  Top  panel:  Schematic  of  the  grating  based 
bistability  scheme. 


Figure  23  shows  the 
wavelength  dependence  of  the 
amplifier  gain  for  different 
modal  gain  powers  and  for 
different  values  of  a .  The 
spectral  shift  can  be  understood 
by  increasing  the  gain  results  in 
a  decrease  in  refractive  index. 

Therefore,  the  transmission 
shifts  to  higher  values. 

Nonlinearities  exist  in  all  gain 
media  as  gain  saturation  occurs 
in  the  large-signal  regime  and  a 
nonlinear  change  in  the  refractive  index 
occurs  along  with  it.  This  nonlinearity, 
together  with,  the  feedback  mechanism,  can 
produce  regions  of  bistability. 

In  order  to  incorporate  these  effects  into  the 
calculations,  an  intensity  dependent  gain  of 
the  form  g  =  g0  /(l  +  ///sat)  is  used.  Since  the 

transfer  matrix  method  uses  g  to  calculate  the 
intensity,  a  self-consistent  approach  is  need. 
The  approach  taken  is  to  use  a  first  order 
approximation  g  =  g0,  and  calculate  the 

intensity  distribution.  Using  the  calculated 
intensity  distribution,  a  new  value  of  g  is 
calculated.  Now  using  this  value  of  g,  the 


81- 

Fig.  23  Gain  dependence  on  wavelength. 


Fig.  24.  Bistability  curves  for  different  gain  values 
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intensity  distribution  is  again  calculated.  This  process  is  repeated  until  the  intensity  distribution 
converges. 

Bistability  exists  for  wavelengths  higher  than  the  transmission  peak.  Hysteresis  curves  are 
plotted  for  values  of  8L  in  the  vicinity  of  the  onset  of  bisability  (8L  =  \  .5 ).  Based  on  realistic 
estimates  of  gain  coefficient  of  colloidal  QDs  and  the  typical  device  length  (500  microns  -  1000 
microns),  we  estimated  the  threshold  power  levels  to  be  in  mW  range  which  can  be  easily 
obtained. 

We  also  carried  out  pump  power  dependent  measurements  on  QDs  spin  coated  on  gratings  that 
were  fabricated  via  soft  lithography.  SEM  images  of  gratings  fabricated  using  electron  beam 
lithography  are  shown  in  Fig.  25.  The  ASE  measurements  were  carried  out  on  gratings  fabricated 


Fig.  25.  Scanning  electron  microscope  image  of  the  gratings  fabricated  on  PMMA  using  electron  beam 
lithography  and  (b)  amplified  spontaneous  emission  from  CdSe  quantum  dots  deposited  on  top  of  a  grating. 
At  high  pump  power  levels,  multimode  ASE  spectrum  is  observed. 


using  softlithography.  Here  we  see  evidence  of  the  different 
modes  arising  due  to  the  feedback  from  the  underlying  grating 
modes. 

The  final  phase  of  the  project  was  dedicated  to  the  fabrication  ogf 
low  cost  DFB  lasers.  This  was  carried  out  by  taking  a  commercial 
DVD  and  stamping  PDMS  onto  it  and  creating  a  negative  mould. 
This  was  later  transferred  onto  Si02  via  dry  etching.  Colloidal 
CdSe  QDs  were  deposited  onto  the  Si02  gratings  via  dip  coating. 
An  atomic  Force  Microscope  image  of  the  grating  is  shown  in 
Fig.  26. 


Fig.  26.  AFM  image  of 
grating  etched  onto  Si02. 


The  samples  were  pumped  using  the  NdYAG  laser  and  the  observed  spectra  is  shown  in  Fig. 
27(a)  along  with  the  power  dependence  in  Fig.  27(b).  There  is  significant  spectral  narrowing  (<  1 
nm).  The  Pout  versus  Pin  graph  also  shows  a  clear  threshold  behavior.  We  are  currently  carrying 
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out  coherence  measurements  to  establish  that  this  is  indeed  a  laser  and  not  amplified  spontaneous 
emission. 


Fig.  27.  (a)  Emission  spectrum  of  the  QDs  coated  on  the  DFB  structure  fabricated  via  soft-lithography  showing 
significant  spectral  narrowing  (<lnm).  (b)  Power  dependence  of  the  emission  showing  clear  threshold  behavior. 
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APPENDIX 


Lasing  equations  and  analysis  of  single  mode  operation 

In  order  to  derive  lasing  equations  for  the  structure  under  consideration  we  expand  the  field  in 
the  two-disk  system  in  terms  of  the  modes  defined  in  Eq.  (20)  -  (23)  and  use  the  orthogonality 
condition,  Eq.  (27)  to  derive  equations  for  slow  amplitudes. 

We  will  describe  the  field  and  the  modes  as  two-dimensional  vectors  defined  in  Eq.  (24),  and 
present  the  modal  expansion  in  the  form 


£}=Z2X» 


l,p 


(33) 


where  first  index  goes  over  non-degenerate  modes,  and  the  second  index  sums  over  degenerate 
ones.  Using  definition  of  the  adjoint  vectors  given  in  Eq.  (25)  and  orthogonality  we  derive 


hP  I 


4  7ik2 


r 

||  dr  tT^)  (r  )i5(1)  (r)  + 

||  dr  (r  )P(2)  (r ) 

V  A 

d2  y 

(34) 


where 


P) 


^  p(l)  ^ 
D(  2) 

\r  J 


(35) 


represent  nonlinear  polarization  due  to  the  gain  medium.  Introducing  slow  changing  amplitude 
approximation,  we  have  instead  of  (34) 


(-k  +  km,i)Ei,P(k) 


Ink 


f 


\ 


||  dr  ^  (r)P(1)  (r)  +  JJ  dr  'F™  (r  )P(2)  (r) 

V  A  D2  ) 


(36) 


We  will  assume  that  there  are  two  lasing  modes,  labeled  by  index  /  =  1,2 ,  characterized  by 
frequencies  Kt  and  spatial  distribution  functions.  The  latter  are  comprised  of  two  degenerate 
modes  corresponding  to  each  value  of  index  / : 

<37> 

p= i 

The  total  field  is  presented  as 

\E)  =  ^\Bl(k))S(k-Ki)  (38) 
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Now  we  follow  Ref.  [4],  For  a  given  lasing  mode  i  the  polarization  term  in  Eq.  (36)  can  be 
presented  in  the  following  form 


\\dr'FS(r)Pm(r)  +  \\dr'FS(r)Pm(r)  +  *£!£,  )E 


where 


where  //f’2)(r)  denotes  the  spatial  profile  of  nonlinear  polarization.  In  the  time-domain  Eq.  (36) 
becomes 


-/  1- 


dk.(Kt)ldE 


dk  dt 


r/f’2){ r)  up  to  the  3rd  order  of  nonlinearity  is  defined  as  follows  [4]: 


tf’2)  =  my^D^Al  1  +  2,S|Er|2  Re[D(^. ) 


+^S|N,'2>|2 D,  (K>  ~ )["(*;  )■ + »'(*<, ) 


where 


A  =  — 


2h2y±yn 


(  k-k  V1  (  kXl 

D(k)=  1  —  i - -  ;DJk)=  1  —  i — 


Aw'/’  is  the  population  imbalance  between  the  two  levels  responsible  for  the  lasing  in  each  disk, 
and  k0 ,  yL  ,  and  yn  represent  frequency  of  the  lasing  transition,  dephasing  and  population 

relaxation  rates  of  the  gain  medium,  respectively  and  d  is  the  dipole  matrix  element  of  the 
transition.  First  we  separate  it  into  linear  and  nonlinear  parts: 

^’2)=<i2)+^32)  (* 


where 
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ni =  2iltrM>D(Kl)A  =  -i^-H,D(K,) 
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2^ZK'2f  Re[DK) 


h=  1 


Respectively,  polarization  matrix  is  also  divided  in  two  parts 
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so  that  Eq.  (41)  can  be  rewritten  as 
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where  we  introduced  integrals 


(50) 
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It  is  convenient  to  present  the  modal  functions  in  the  following  form: 

F®  =  R~xeim(pY^\r)\  F™  =  R~l e~im<pYx{1)  (r); 

F®  =  (-1  )mR~le~im,pY^\r);  F$  =  (~\)m  R-xeim<pYf\r) 
F2(j}  =  eim,p R-'Y^  (r);  F2(2)  =  frVmfT2(2)(r); 

F2a2  =  (-\)m  R-'e-,m,pYlu(r);  =  (-1  )m R~leim^2\r) 


where  functions  Y^{ 2)  are  defined  by  comparing  Eq.  (5 1)  to  Eq.  (20)  -  (23).  Most  of  the 

coefficients  turn  zero  because  of  the  integration  over  angle  cp .  The  only  non-zero 

coefficients  are  (Self-satuiration) 
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where 
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For  cross-saturation  terms  we  have  similalry 
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The  system  of  lasing  equations  for  modal  amplitudes  can  now  be  presented  in  the  following  form 
JE, 
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and 


b9  =  to$X  +  tl%'b?  («) 

We  also  neglected  the  derivative  of  the  eigenfrequency  with  respect  to  the  spectral  parameters, 
which  is  allowed  for  high-Q  resonators.  Eq.(56)  and  Eq.  (57)  describe  dynamics  of  the  field  in 
the  first  pair  of  degenerate  modes,  while  Eq.  (58)  and  (59)  deal  with  the  second  pair.  In  the 
case  of  the  single  mode  lasing,  only  one  of  these  pairs  have  non-zero  field.  For  instance,  we  can 
choose  E2  j  =  E2  2  =  0 ,  in  which  case  Eq.  (56)  and  (57)  reduce  to 
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Eq.  (62)  shows  that  the  degenerate  modes  cross-saturate  each  other  with  the  coefficient  of  cross 
saturation  two  times  larger  than  that  of  self-saturation.  Therefore,  assuming  that  one  of  the 
modes  starts  lasing  first  it  will  saturate  population  inversion  for  the  second  mode  suppressing  its 
lasing.  Since  both  degenerate  modes  have  the  same  threshold,  one  can  assume  that  only  one  of 
those  modes  will  play  a  role  in  the  devices  under  construction. 

Therefore,  we  can  assume  that  the  two-mode  lasing  is  only  possible  with  non-degenerate  modes, 
i.e.  Ej  j  and  E2 1 ,  we  can  study  the  mode  dynamics  with  reduced  equations 
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Presenting  field  in  the  form EiX  =  |E. | exp ( -iy/i )  and  separating  real  and  imaginary  parts  of  (63) 
we  have  equations  for  intensity 


and  phase 
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nd 


42=a 


1,2 


a  _jj  2/rk()d2An(]  K^k, 


fin] 

z  ^-nK\d  An0 


(*,-*»)  +y] 

Yl 


tin]  ( Ki-k0 )  +rl 


(64) 


(65) 


(66) 


(67) 


The  effective  coupling  coefficient  determining  stability  of  various  lasing  regimes  in  this  case  is 
given  by 


q  _  Re[42^i2]Re[^2Ai] 
Re  [  \  \bx  j  ]  Re  [  h22b22  ] 


(68) 
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Neglecting  population  pulsation  terms  in  (60)  we  have 

BpC,)Rep>(^)] 


and 


hn 


h22  — 


\i 


2  7ik0 

1  ( 

nd 

2  7rk0 

1 

n] 

2  7ik0 

1  | 

n] 

hY\\ 

f  d2  ^ 


yhyLRj 
f 


\2 


h2l  —  ■ 


irk,,  1 


KhyLRj 

Khy±Rj 


(  J  2  \ 


"i  »r, 


!>r  R  > 


D(^)Re[D(^,)] 

D(^,)Re[D(A:,)] 


C  = 


Re  \_D(K{)  bn  ]  Re[Z)  (K2  )b2l] 
Re[D(^)6ll]Re[D(A'2)A2] 


(69) 


(70) 


Behavior  of  the  laser  depends  on  the  parameters  g,  ( Kx ),  ( K2 ) :  the  lasing  begins,  when  one  of 

these  parameters  becomes  greater  than  zero.  In  the  single  mode  regime,  when  only  one  mode 
lazes,  we  have 

t=V0,(K,>+R#AA) 

Consider 


b  = 


2 irk,  1  (  d2 
nd  hY\  f'Y_R  - 


D(A:,)Re[D(A:,)]; 


Re[D(Jf,)] 

Im[£>(4)] 


i 


Yl 


(^i)2 


OO  +(*,-*»)' 


=  g(^) 


002+(*,-*o)2 


=/(*,) 


(72) 


Lasing  threshold  is  defined  as 

An(0  =  (^,  -^o)2+ri  hrfrl  _  k„ 

0  Ink^d2  2nk0d2 

where  we  used  equation  for  the  lasing  frequency  at  the  threshold 
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r,+K(K,) 

,  rMK,)+KK(^)  , 

”  rL+V(K,)  "  r^k’(K,) 


(74) 


It  is  important  to  determine,  which  mode  has  a  lower  threshold,  while  the  absolute  threshold 
value  is  not  important.  Thus  we  can  define  the  “threshold  coefficient” 


_  Hn 

Ink^d2  R 


Defining  thresholds  normalized  by  y/th  we  have 


=  Yi 


1  + 


(*<~*o)2 

(r±+Yi)2 


(75) 


(76) 


Where  all  xi  and  y.  are  dimensionless  real  and  imaginary  parts  of  the  eigenfrequency  of  the 
respective  mode,  calculated  at  the  lasing  frequency,  x0  and  yL  are  dimensionless  atomic 

frequency  and  gain  width  respectively.  Dimensionless  equation  for  the  threshold  lasing 
frequency  is 


_yixt{Xi)  +  xji(Xi) 

r±+/i(Xi ) 

In  the  stationary  single  mode  state: 

/  few  Ryyrjt  a»,-a<»> 

'  Re[»A]  d2  Re[>„Z>(tf,)] 

-  An, -A  ,  d 2 

'  Re[6,D(X,)]’  ‘R2t?nr, 


(77) 


(78) 


where  we  introduced  dimensionless  intensity.  For  the  stationary  single  mode  frequency  we  have 
equation: 


Re  [6,4] 


Im 

\D(x,y_ 

Re 

\D(X,)_ 

where  cr  and  ^  can  be  rewritten  as 


ai=k'i-An0f(Xi) 

=(a«0  -A4)g(W) 


(79) 


(80) 
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and 


D(k)=g(x,)+nx, ) 


(/,)  +(X,-x o)' 


-;/(*,)  = 


(rJ2+(^-^o)2 


(81) 


The  derived  formulas  lay  foundation  for  the  analysis  of  laser’s  behavior  for  various  parameters 
of  the  system  and  for  design  of  the  desirable  bistable  behavior. 
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